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Abstract 

Wc give an ©(nlog'^n) algorithm that, given an n-nodo directed planar graph with arc capacities, 
a set of source nodes, and a set of sink nodes, finds a maximum flow from the sources to the sinks. 
Previously, the fastest algorithms known for this problem were those for general graphs. 



1 Introduction 



The maximum flow problem with multiple sources and sinks in a directed graph with arc-capacities is, 
informally, to find a way to route a single commodity from a given set of sources to a given set of sinks 
such that the total amount of the commodity that is delivered to the sinks is maximum subject to each arc 
carrying no more than its capacity. In this paper we study this problem in planar graphs. 

The study of maximum flow in planar graphs has a long history. In 1956, Ford and Fulkerson introduced 
the max st-flow problem, gave a generic augmenting-path algorithm, and also gave a particular augmenting- 
path algorithm for the case of a planar graph where s and t are on the same face Researchers have since 
published many algorithmic results proving running-time bounds on max st-flow for (a) planar graphs where 
s and t are on the same face, (b) undirected planar graphs where s and t are arbitrary, and (c) directed 
planar graphs where s and t are arbitrary. The best bounds known are (a) 0{n) 14 , (b) O(nloglogn) |17) . 
and (c) O(nlogn) [5^, where n is the number of nodes in the graph. 

Maximum flow in planar graphs with multiple sources and sinks was studied by Miller and Naor |29) . 
When it is known how much of the commodity is produced/consumed at each source and each sink, finding 
a consistent routing of flow that respects arc capacities can be reduced to negative-length shortest paths, 
which we now know can be solved in planar graphs in 0(nlog^ n/ log log n) time [301. Otherwise, Miller and 
Naor gave an 0{n\og'^^^ n) algorithm for the case where all the sinks and the sources are on the boundary 
of a single face, and generalized it to an 0{k^n^^'^ log^ n)-time algorithm for the case where the sources and 
the sinks reside on the boundaries of k different facesQ 

However, the problem of maximum flow with multiple sources and sinks in planar graphs without any 
additional restrictions remained open. In general (i.e., non-planar) graphs, multiple sources and sinks can 
be reduced to the single-source single-sink case by introducing an artificial source and sink and connecting 
them to all the sources and sinks, respectively — but this reduction does not preserve planarity. For more 
than twenty years since the problem was explicitly stated and considered |29| . the fastest known algorithm 
for computing multiple-source multiple-sink max-flow in a planar graph has been to use this reduction in 
conjunction with a general maximum-flow algorithm such as that of Sleator and Tarjan |32) which leads 
to a running time of 0(n^ log n). For integer capacities less than U, one could instead use the algorithm 
of Goldberg and Rao [TU], which leads to a running time of 0(n^-^ log n log C/). No planarity-exploiting 
algorithm was known for the problem. 

Our Result The main result of this paper is an algorithm for the problem that is optimal up to a small 
poly-logarithmic factor. 

Theorem 1.1. There exists an algorithm that solves the maximum flow problem with multiple sources and 
sinks in an n-node directed planar graph in 0(7ilog'^ n) time. 

Application to computer vision problems Multiple-source multiple-sink min-cut arises in addressing 
a family of problems associated with the terms metric labeling (Kleinberg and Tardos, H^), Markov Ran- 
dom Fields H], and Potts Model (see also [HITS])- In low- level vision problems such as image restoration, 
segmentation, stereo, and motion, the goal is to assign labels from a set to pixels so as to minimize a penalty 
function. The penalty function is a sum of two parts. One part, the data component, has a term for each 
pixel; the cost depends on the discrepancy between the observed data for the pixel and the label chosen for 
it. The other part, the smoothing component, penalizes neighboring pixels that are assigned different labels. 

For the binary case (when the set of available labels has size two), finding the optimal solution is reducible 
to multiple-source multiple-sink min-cut. [12j . For the case of more than two labels, there is a powerful and 
effective heuristic [4] using very-large-neighborhood [1] local search; the inner loop consists of solving the 
two-label case. The running time for solving the two-label case is therefore quite important. For this reason, 
researchers in computer vision have proposed new algorithms for max flow and done experimental studies 
comparing the run-times of different max-flow algorithms on the instances arising in this context [31 [27] . All 
of this is evidence for the importance of the problem. 

^The time bound of the first algorithm can be improved to 0(n log n) using the linear-time shortest-path algorithm of 
Henzinger et al. |14) . and the time bound of the second algorithm can be improved to 0(fc^?ilog'^ n) using the 0(ra log n)-time 
single-source single-sink maximum flow algorithm of Borradaile and Klein j2,. 
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For the (common) case where the underlying graph of pixels is the two-dimensional grid graph, our result 
yields a theoretical speed-up for this important computer-vision subroutine]^ 

Hochbaum |15j describes a special case of the penalty function in which the data component is convex 
and the smoothing component is linear; in this case, she shows that an optimal solution can be found in 
time 0{T{m,n) + nlogU) where U is the maximum label, and T(m,n) is the time for finding a minimum 
cut. She mentions specifically image segmentation, for which the graph is planar. For this case, by using 
our algorithm, the optimal solution can be found in nearly linear time 



Application to maximum bipartite matching Consider the problem of maximum matching in a bipar- 
tite planar graph. It is well-known how to reduce this problem to multiple-source, multiple-sink maximum 
flow. Our result is the first planarity-exploiting algorithm for this problem (and the first near- linear one). 



Techniques To obtain our result, we employ a wide range of sophisticated algorithmic techniques for 
planar graphs, some of which we adapted to our needs while others are used unchanged. Our algorithm 
uses pseudoflows [TTl [TS] and a divide-and-conquer scheme influenced by that of [TH] and that of [2^] . We 
adapt a method for using shortest paths to solve max st-flow when s and t are adjacent [13], and a data 
structure for implementing Dijkstra in a dense distance graph derived from a planar graph [8]. Among the 
other techniques we employ are: using cycle separators )28j recursively while keeping the boundary nodes on 
a constant number of faces [521 [SSI [H] , an algorithm for single-source single-sink max flow [H [3 1 an algorithm 
for computing multiple-source shortest paths [231 [S], a method for cancelling cycles of flow in a planar 
graph |20J , an algorithm for finding shortest paths in planar directed graphs with negative lengths [23] [30] , 
and a data structure for range queries in a Monge matrix [21", . 



1.1 Preliminaries 

We assume the reader is familiar with the basic definitions of planar embedded graphs and their duals (cf. 
[3j). Let G = {V,E) be a planar embedded graph with node-set V and arc-set E. We use the term arc to 
emphasize that edges are directed. The term edge is used when the direction of an arc is not important. For 
each arc a in the arc-set E, we define two oppositely directed darts, one in the same orientation as a (which 
we sometimes identify with a) and one in the opposite orientation [2 . We define rev(-) to be the function 
that takes each dart to the corresponding dart in the opposite direction. It is notationally convenient to 
equate the edges, arcs and darts of G with the edges, arcs and darts of the dual G* . 

Let S GV he a. set of nodes called sources, and let T C y — 5' be a set of nodes called sinks. 

A flow assignment /(•) is a real- valued function on darts satisfying antisymmetry: 

firevid)) = -f{d) 
A capacity assignment c(-) is a real- valued function on darts. 

A flow assignment /(•) respects the capacity of dart d if f{d) < c{d). /(•) is called a pseudoflow if it 
respects the capacities of all darts. 

For a given flow assignment /(•), the net inflow (or just inflow) node v is inflow j^(ti) = X^dart d-hcad{d)=v /("^Q 
The outflow of v is outflow^(w) = — inflow/(w). A flow assignment /(•) is said to obey conservation at node 
V if inflow j^(ti) = 0. A feasible flow is a pseudoflow that obeys conservation at every node other than the 
sources and sinks. A feasible circulation is a pseudoflow that obeys conservation at all nodes. The value of 
a feasible flow /(•) is the sum of inflow at the sinks, X^tGT inflow ^'(i) or, equivalently, the sum of outflow at 
the sources. The maximum flow problem is that of flnding a feasible flow with maximum value. 

For two flow assignments /, /', the addition f + f is the flow that assigns f{d) + f'{d) to every dart d. 

A residual path in G is a path whose darts all have strictly positive capacities. For two sets of nodes 
A, B, A ^ B is used to denote the existence of some residual a-to-6 path in G for some nodes a £ A and 
b £ B. Conversely, A ^ B is used to denote that no such path exists. We will omit the graph G when it is 
clear from the context. 



^Note that the single-source, single-sink max-fiow algorithm of ^2^ was implemented by computer-vision researchers | 31| and 
found to be useful in computer vision and to be faster than the competitors. 

3An equivalent definition, in terms of arcs, is infiow_f (j;) = EaeE:head(a)=i, ' Eae£;:taii(a)=« /(«)■ 
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The residual graph of G with respect to a flow assignment /(•) is the graph Gf with the same arc-set, node- 
set, sources and sinks, and with capacity assignment Cf{-) such that for every dart d, Cf{d) = c{d) — f{d). 

It is weh known that a feasible flow / in G is maximum if and only if S* ^ T. 

Let / be a pseudoflow in a planar graph G. Let denote the set of nodes {v G V — (S Li T) : 
inflow/(t') > 0}. Similarly, let denote the set of nodes {v E V — {S U T) : inflow/(i') < 0}. Suppose 

S U V'^ S> T U V~ . For a graph with n nodes and m edges, there exists an 0(r7ilogn)-time algorithm that 
converts the pseudoflow / into a maximum feasible flow /' [111132]. In planar graphs, this can be done in 
linear time by first canceling flow cycles using the technique of Kaplan and Nussbaum [20 , and then by 
sending back flow from and into V~ in topological sort order. See Appendix [c| for details. 



1.2 Overview of the algorithm 

Consider the following recursive approach for finding a maximum multiple-source multiple-sink flow: split 
the input graph G in two using a simple cycle separator G |28| and recursively solve the max flow problem 
in the two subgraphs. When the two recursive calls have been executed, in each of the two subgraphs there 
is no residual path from any source to any sink. If we further make sure that in each of the two subgraphs 
there is no residual path from any source to C and from G to any sink, then, since C is a separator, there 
is no residual path from any source to any sink in G. 

We therefore solve a slightly more general problem recursively in the two subgraphs: roughly speaking, 
find a fiow such that there is no residual path from a source to a sink or to G and no residual path from G 
to a sink (Section [2]). After the two recursive calls there is no residual path from any source to any sink in 
G. However, the requirement that there is no residual path from any source to G and from G to any sink 
cannot be achieved by a feasible flow but rather by a pseudoflow in which there might be excess inflow or 
excess outflow on nodes of G. We deal with this by solving a new max flow problem where nodes of G are 
treated as sources and sinks, limited in supply/demand by their excess inflow/outflow (Section [3|. 

We exploit a relation between primal circulations and dual shortest paths to maintain a succinct repre- 
sentation of the flow during critical phases of the algorithm, using the fact that there are only 0{y/n) sources 



and sinks, all cyclically ordered on G (Section 3.2). Even though our representation does not explicitly store 
the flow on nearly any arc in the graph, we can augment it efficiently towards optimality while maintaining 
feasibility. An important tool we use is Fakcharoenphol and Rao's efficient implementation of Dijkstra's 
algorithm [8 , which we adapt to our needs. 

The resulting pseudoflow can then be turned into a max flow in linear time using existing techniques 
(Appendix O). This leads to an O(nlog^n) time algorithm for max flow. 



2 The Algorithm 

Now we describe the algorithm, referred to as MultipleSourceMultipleSinkMaxFlow, in more detail. 
In order to treat nodes of the cycle separator both as sources and as sinks in recursive calls, we introduce a 
new node set A. At the top recursion level, A — (d. In general, A has constant size; more precisely \A\ < 6. 

MultipleSourceMultipleSinkMaxFlow(G, c,S,T,A) 

Input: a directed planar graph G with non-negative capacities c, a set S of source nodes, a set T of sink 
nodes, a set A of at most six nodes 

Output: a pseudoflow / obeying conservation everywhere but S, T, A and s.t. S^T,S^A, A^T. 

The algorithm finds a simple cycle separator G pS] and contracts all edges of G except one. This, 
essentially, merges all the nodes of G into a single supernode v, and turns G into a self loop. For simplicity 
of presentation we assume that no sources or sinks lie on The algorithm then recursively solves the 
problem on the subgraphs enclosed and not enclosed by that self loop (the self loop itself need not be 
included in any of the subgraphs), adding the supernode v that represents G to the set A. In order to keep 

^This does not lose generality since if a node u £ C is a source, one can introduce a new node u' and an arc u'u whose 
capacity equals the sum of capacities of arcs outgoing from u. Now consider n' as a source instead of u. Since the separator 
has just O(v'n) nodes, this will not affect the running time. Sinks can be handled in a similar fashion. 
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the cardinality of A at most six, the algorithm alternately applies the cycle separator theorem with weights 
uniformly distributed on all nodes and uniformly distributed on only the nodes of A. This technique is 
similar to that used in [251 Ell H] • 

After the recursive calls, the algorithm uncontracts the edges of C. At this stage there are no residual 
paths between sources and sinks in the entire graph, but there might be excess inflow (positive or negative) 
at the nodes of C. The algorithm then calls the procedure FixConservationOnPath that pushes flow 
between the nodes of C so that there are no residual paths between nodes of C with positive inflow and 
nodes of C with negative inflow (the path in the name of the procedure is the cycle C without one edge). 
This procedure is discussed in Section |3l the interface is: 

FixConservationOnPath(G', P, c, /g) 

Input: a directed planar graph G, simple path P, capacity function c, and a pseudoflow /o 
Output: a pseudoflow / s.t. (i) f ~ fo satisfies conservation everywhere but P, and 

(a) {v & P : inflow('i;) > 0} ^ {w G P : inflow(u) < 0}. 
Running Time: 0(n log^ n/ loglog n + jPplog^n) 

Next, the algorithm iterates over the nodes of A. The algorithm calls the procedure CycleToSin- 
GLeSinkLimitedMaxFlow that, roughly speaking, pushes as much excess flow as possible from C to a^. 
If Cf is the set of nodes of C that are reachable via residual paths from some node of C with positive inflow 
at the beginning of iteration z, CycleToSingleSinkLimitedMaxFlow pushes flow among the nodes of 
Cf and from the nodes of to ai. The result is that remaining inflow at nodes of Cf is non- negative and 
there are no residual paths from nodes of C with positive inflow to at. See Section |4] the interface is: 

CyCLEToSiNGLESiNKLIMITEDMAXFLOW(G, C, /q, G, flj) 

Input: a directed planar graph G with capacities c, a pseudoflow /o, a simple cycle G, a sink a^. 

Assumes: Vw e G+, infiow/Q(t;) > 0, where G+ = {w G G : {a; G G : inflow j^^, (a;) > 0} -4' v}. 

Output: a pseudoflow / s.t. (i) f — fo obeys conservation everywhere but G+ U {t}, (ii) Vv E 

Gf 

G+, inflow^(w) > 0, (in) {v e C : inflow ^^(w) > 0} ^ a^. 
Running Time: 0(n log^ n/ log log n + jGplog^n). 

A similar procedure SingleSourceToCycleLimitedMaxFlow is called to push flow from to G to 
eliminate as much negative inflow as possible (using a similarly defined set C^). 

Finally, the algorithm pushes back flow from any nodes of G with positive inflow to S and pushes flow 
back from T into any nodes of G with negative inflow. 

Algorithm 1 MultipleSourceMultipleSinkMaxFlow(G, c,S,T,A) 

Input: a directed planar graph G with non-negative capacities c, a set S of source nodes, a set T of sink 
nodes, a set A consisting of a constant number of nodes A = {ai}i^i- 

G G G 

Output: a pseudoflow / obeying conservation everywhere except S,T,A and s.t. S-^^TjS-y^A, A-y^T. 

1: add zero capacity arcs to triangulate and 2-connect G (required for simple cycle separators) 
2: find a balanced (w.r.t. |G| and \A\, alternately) cycle separator G in G disjoint from S and T 
3: let P be a path comprising of all of G's edges except one edge e 

4; contract all the edges of P, turning e into a self loop incident to the only remaining node u of G 
5; let Gi and G2 be the subgraph of G enclosed and not enclosed by e, respectively 

6; / := MULTIPLESOURCEMULTIPLESlNKMAXFLOw(Gi,C|Gi,S'n Gi,rn Gi, (An Gi) U {u}) 

7; / := / + MULTIPLESOURCEMULTIPLESlNKMAXFLOW(G2,C|G2,S'nG2,PnG2,(AnG2)U{t;}) 

8; uncontract the edges of P 

9; / := FixConservationOnPath(G, P, c, /) 

10; for i = 1, 2, . . . , fc 

11; / := CycleToSingleSinkLimitedMaxFlow(G, c, /, G, at) 

12; / := SiNGLESOURCEToCyCLeLiMITEDMAxFLOW(G, C, /, Qi, C) 

13; push positive excess from G to S* and negative excess to G from T 
14; return / 
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Correctness of Algorithm [T] The correctness of Algorithm [T] is proved in detail in Appendix |Xj The 
proof consists of a sequence of lemmas that prove that each step of the algorithm eliminates some undesired 
residual paths without reintroducing undesired residual paths. The arguments used are elementary. 

Running Time Analysis The number of nodes of the separator cycle C used to partition G into Gi 
and G2 is 0{^J\G\). Therefore, each invocation of FixConservationOnPath, CycleToSingleSinkLim- 
itedMaxFlow and SingleSourceToCycleLimitedMaxFlow in G takes 0(|G| log^ |G|/loglog |G| + 
|Gplog2|G|)=0(|G|log^|G|)time. 

The way we recursively partition into subgraphs is very similar to that of Fakcharoenphol and Rao ^ . 
In their algorithm, they spend 0(|G'| log^ |G'|) time on each subgraph G' in the recursive decomposition and 
prove a total time bound of 0{n\o^ n). By the same arguments, our algorithm runs in 0{n\og^ n) time. 

3 Eliminating Residual Paths Between Nodes with Positive Inflow 
and Nodes with Negative Inflow on a Path 

In this section we present an efhcient implementation of the fixing procedure which, roughly speaking, given 
a path with nodes having positive, negative, or zero inflow, pushes flow between the nodes of the path so 
that eventually there are no residual paths from nodes with positive inflow to nodes with negative inflow. 

We begin by describing an abstract algorithm for the flxing procedure. The abstract algorithm is given 
as Algorithm [2] It is similar to a technique used by Venkatesan and Johnson [TS] . Let M be the sum of 
capacities of all of the darts of G. The algorithm first increases the capacities of darts of the path P and 
their reverses by M. Let pi,p2, ■ ■ ■ ,Pk+i be the nodes of P. The algorithm processes the nodes of P one after 
the other. Processing pi consists of decreasing the capacities of di — PiPi+i and rev{di) by M (i.e., back to 
their original capacities) , and trying to eliminate positive inflow x at pi by pushing at most x units of flow 
from Pi to Pi+i- The intuition for doing so is that the flow after the push either obeys conservation at pi or 
there are no residual paths from pi to any of the other nodes of P (this is where we use the large capacities 
on the darts between unprocessed nodes). See appendix |b] for a formal proof of correctness. Negative inflow 
at Pi is handled in a similar manner by pushing flow from Pi+i to pi. 



Algorithm 2 AbstractFixConservationOnPath(G, P, c, /q) 

Input: directed planar graph G, simple path P = did2 ■ ■ ■ dfc, capacity function c, and pseudoflow /o 
Output: a pseudoflow /' s.t. (i) f — fo satisfies conservation at nodes not on P, and (ii) with respect to 
/', there are no residual paths from nodes of P with positive inflow to nodes of P with negative inflow. 

1: /' - fo 

2: c[d] = c[d] + M for all darts d of P U rev(P) 

3: for i — 1,2, . . . ,k 

4: let Pi and Pi+i be the tail and head of di, respectively 

% reduce the capacities of d and rev(d) by M and adjust the flow appropriately 
5: for d e {di,rev{di)} 
6: c[d] := c[d] - M 
7: f[d]:=mm{f'[d],c[d]} 
8: /'[rev(d)] := -f[d] 

9: excess := inflow at pi 

10: if excess > then d := di else d :— rev(di) find in which direction flow should be pushed 

11; add to /' a maximum feasible flow from tail(c?) to head((i) with limit excess 
12: return /' 
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3.1 An InefRcient Implementation 

In this section, we give an inefficient implementation of line [Tl] of the abstract algorithm. This will facilitate 
the explanation of the efficient procedure in the next section. We first review the necessary ideas and tools. 

3.1.1 Hassin's algorithm for maximum st-planar flow 

An st-planar graph is a planar graph in which nodes s and t are incident to the same face. Hassin [13^ gave 
an algorithm for computing a maximum flow from s to t in an st-planar graph. We briefly describe this 
algorithm here since we use it in implementing line [TT] of Algorithm [2] 

Hassin's algorithm starts by adding to G an artificial infinite capacity arc a from t to s. Let d be the 
dart that corresponds to a and whose head is t. Let t* be the head in G* of the dual of d. Compute in the 
dual G* a shortest path tree rooted at t* , where the length of a dual dart is defined as the capacity of the 
primal dart. Let denote the shortest path distances from t* in G* . Consider the flow 

p[d'] = </)[headG. {d')] - ^[tailc- {d')]foT all darts d' (1) 

After removing the artificial arc a from G, p is a maximum feasible flow from s to i in G. We say that (p is 
a face potential vector that induces p. 

In our algorithm we are interested in a max flow with limit x from s to t rather than a maximum flow, 
i.e., a flow whose value is at most a given number x but is otherwise maximal. It is not difficult to see that 
setting the capacity of the artificial arc a to a; instead of infinity results in the desired limited max flow [22] . 

In our implementation, instead of using an artificial arc from t to s, we use an existing arc whose endpoints 
are s and t as the arc a above. In order for this to work, the capacity of the dart d that corresponds to a 
and whose head is t must be zero (as is indeed the case if a is an artificial arc from t to s). This can always 
be achieved by first pushing flow on d to saturate it. Also note that in this case, we do not remove a from 
G. Hence, p is a feasible circulation, rather than a maximum flow, since flow is being pushed back from t to 
s along a. To convert it into a maximum flow one just has to undo the flow on a. If we define / by 

-r _ J ~p[d'] if d' corresponds to a , , 

" \ otherwise ' ^ ' 

then p + / is a maximum feasible st-flow. We will use the fact that this flow can be represented implicitly 
by the face potential vector ^, and the flow values f[d'] for the two darts corresponding to a. 

3.1.2 The Inefficient Implementation 

Recall that /o is the flow at the beginning of the procedure. Observe that the change to the flow in iteration i 



of the abstract algorithm (line II ) is a flow between the endpoints of the dart di. As discussed in Section 3.1.1 
this flow can be represented as the sum of [i) a circulation pi and {ii) a flow on dj and rev((ii). Summing 
over the first i iterations, the flow /' at that time can be represented as the sum 

f' = P + f (3) 

where p = Pj ^ circulation and / is a flow assignment that differs from /o only on the darts of 

{dj}j=i and their reverses. 

The inefficient implementation of line [TT] of the abstract algorithm appears as Algorithm [3] We now 
describe it. The total flow /' is maintained by representing / and the circulation p as in Eq. ([3|. / is 
represented explicitly, but, in preparation for the efhcient implementation, the circulations Pj are represented 
implicitly by the face-potentials cj)j. By linearity of Eq. ([!]), the sum (p — X]j=i 't'j is the face potential vector 
that induces the circulation p. 



Recall that d is the dart of C such that flow needs to be sent from tail(d) to head(d) (line 10 of Algo- 
rithm [2]) . In lines [l] ~ |3j the procedure pushes as much as possible on d itself. Consequently, either d is 
saturated or conservation at pi is achieved. 

Next, an implementation of Hassin's algorithm pushes a maximum flow with limit |inflow(pi)| from tail((i) 
to head{d). The procedure first sets the length of darts to their residual capacities (lineji]) and sets the length 
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of rev((i) to be the flow limit |inflow(pi)| (linejs]). Since the flow maintained is feasible, all residual capacities 
are non-negative. The procedure then computes all the headc* (d)-to-/ distances cf)^[f] in G* using Dijkstra's 
algorithm (line|6]). Let pi be the circulation corresponding to the face-potential vector cj}^. The procedure 
sets to/ equal to p[d\ in line]?] then subtracts val bom f[d\ and adds it to /[rev(d)]. Finally, in the last line, 
the current circulation is added to the accumulated circulation by adding the potential 0j to (p. 



Algorithm 3 Incflicicnt Implementation of line 11 of AbstractFixConservationOnPath (Algorithm 2) 



% push flow on d to make its residual capacity zero as required for Hassin's algorithm 
residual capacity :— c[d] — f[d] — (p[headG* (d)] + 0[tailG*(rf)] 

val := mm{residual capacity, |inflow(pi)|} amount of flow to push on d 

f[d] := f[d] + val ; /[rev(d)] := -f[d] 

% push excess inflow from tail(d) to head(d) using Hassin's algorithm 

let £[d'] :— c[d'] — f[d'] — 0[headG* (rf')] + 4>[^^^^G*{d')] for all darts d' E G lengths arc residual capacities 

^[rev(d)] |inflow(p.i)| set the limit on the residual capacity of Tcv{d) 

face potential are distances in G* from hcad^* {d) w.r.t. residual capacities 
the amount of flow assigned to d by the circulation corresponding to 0^ 
f[d] :— f[d] — val ; /[rev((i)] := — /[(d)] do not push the circulation on d and rev(ci) 

(p ^ (p -\- (p^ accumulate the current circulat 



ion 



3.2 An Efficient Implementation 

In this section we give an efficient implementation of Algorithm [2j We first review the necessary tools. 

3.2.1 Fakcharoenphol and Rao's Efficient Dijkstra Implementation 

Let A be a set of nodes. Let H he a. planar graph in which the nodes of X lie a single face. Let xi,X2, - ■ ■ 
be the clockwise order of the nodes of X on that face. Let P be a set of darts not necessarily in the graph H 
whose endpoints are nodes in X. Fakcharoenphol and Rao [8] described a data structure that can be used 
in a procedure that efficiently implements Dijkstra's algorithm in H U P. The procedure takes as input a 
table D such that D[i, j] stores the distance between Xi and Xj in H, an array £ that stores the lengths of the 
darts in P, and a node w e A. It is assumed that the lengths in D and in £ are non-negative. The procedure 
outputs the distances of the nodes of A from v in H U P in 0{\X\ log^ |A| -I- \P\ log |A|)-time. 

We mention a technical issue whose importance will become apparent in the sequel. The procedure 
partitions the table D into several subtables {Dq,}q that correspond to distances between pairs of disjoint 
sets of nodes of A, where each set consists of nodes that are consecutive in A. It is assumed in [8] footnote 
on p. 884] that for each such subtable Dq, a data structure that supports range minimum queries of the form 
niinj-^<cj<j2{Da[i, j]} for every i,ji,j2 is given. Fakcharoenphol and Rao note that such a data structure can 
be easily implemented by using a range-search tree [6 for every row i of Dq . The time required to construct 
all of the range-search trees for Dq, is proportional to the size of D^. 

3.2.2 Price Functions, Reduced Lengths and FR-Dijkstra 

For a directed graph G with dart-lengths £{■), a price function is a function from the nodes of G to the 
reals. For a dart d, the reduced length with respect to p is £^{d) — £{d) + p{ta,i\{d)) — p{heeLd(d)). A feasible 
price function is a price function that induces nonnegative reduced lengths on all darts of G (see |18|). 

Single-source distances form a feasible price function. Suppose that, for some node r of G, for every node 
V of G, p{v) is the r-io-v distance in G with respect to £{■). Then for every arc uv, p{v) < p{u) + £{uv), so 
^(t){uv) > 0. Here we assume, without loss of generality, that all distances are finite (i.e., that all nodes are 
reachable from r) since we can always add arcs with sufficiently large lengths to make all nodes reachable 
without affecting the shortest paths in the graph. 

We will use the following variant of Fakcharoenphol and Rao's efficient Dijkstra implementation. The 
procedure FR(D,f, cp"^ , v) takes as input the table D and the array £ as described above, as well as a feasible 
price function (p'^ on the nodes of A and a node v G X. It outputs the distances of the nodes of A from v 
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Figure 1: An example illustrating that the nodes of X* are 
on the boundary of a single face of H* . The diagram shows 
part of the graph G and some edges of its dual G* . Edges of 
G are solid blue. Edges of P are bold. Dual edges are double 
lined red. The dual edges of P are in double lined dashed 
red. 



in H U P w.r.t. the reduced lengths w.r.t. (p'^ . We stress that lengths in D and in £ may be negative, but 
the reduced lengths are all non-negative. The computa tion ta kes 0(|X| log^ \X\ + |P| log \X\) time. 



This differs from the procedure described in Section 3.2.1 only in the existence of the price function (p'^ . 



We cannot afford to compute the entire table of reduced lengths since that would dominate the running 
time of the algorithm in [8 . Instead, whenever their algorithm requires some specific reduced length, we 
can compute it in constant time from D. This, however, does not suffice. Recall that the algorithm in [5] 
requires that, for each of the subtables D^, range-search trees that support range minimum queries of the 
form minjj <j<j2 {Dq [i, j] + [xi] — (f>^ [xj\ } for every z, ji , j2 are given. Note that the results of such queries 
may be different for different price functions. Computing the range-search trees would take 0(|Xp) which 
will dominate the running time of the entire procedure. To overcome this difficulty we use Monge range-query 
data structures, due to Kaplan and Sharir [21 , which can be constructed from the table D in 0(|^| log \X\) 
time, and answer queries of the desired form in 0(log \X\) time. 

3.2.3 The Procedure 

Finally, we describe the efficient implementation. The procedure keeps track of the infiow at each node of P 
in an array v[-]. As in the inefficient implementation, the procedure will maintain the total flow as the sum 
of a circulation p and a flow assignment / that differs from /o only on the darts of P U rev(P). Initially / 
is set equal to /q. The circulation p will be represented by a face-potential vector 0. However, we will show 
that it suffices to maintain just those entries of that correspond to faces incident to P. 

Define each dart d's length by l{d) = c[d] — f[d]. Let X* be the set of endpoints in the planar dual G* 
of the darts of P (i.e. the primal faces incident to P). Let H* be the graph obtained from G* by removing 
the darts of P. Note that in H* , all the nodes of X* that did not disappear (i.e., that have degree greater 
than zero) are on the boundary of a single face, see Figure [l] 

The procedure uses a multiple-source shortest paths (MSSP) algorithm [231 [5] to compute a table D*[-, •] 
of X*-io-X* distances in H* with respect to the lengths £{■). 

The implementation of lines [I]- [lO] of the abstract algorithm using the chosen representation of /' is 
straightforward. We therefore focus on the implementation of line IT of Algorithm [2] given as Algorithm |4] 



Consider iteration i of the algorithm. The main difference between the inefficient implementation and 
the efficient one is in implementing the Dijkstra step for computing shortest paths in the dual. Instead of 
computing the entire shortest-path tree, the procedure computes just the distance labels of nodes in X* . 
This is done using the FR data structure, whose initialization requires the X*-to-X* distances in H* with 
respect to the residual capacities. We now explain how these distances can be obtained. 
The flow on a dart d in the primal is 

f[d] + 0[headG. (d)] - 0[tailG* {d)]. (4) 

Therefore the residual capacity of d is 

{c[d] - f[d]) - 0[headG* {d)] + 0[tailG. {d)] (5) 
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which is its reduced length £(j,[d] with respect to the length £(■) and price function </>. 

Suppose that d belongs to H* , i.e. d is not one of the darts of PUrev(P). The procedure never changes 
f[d], so f[d] — fo[d]. Therefore £{d) — c[d] — fold]. These lengths are known at the beginning of the 
procedure's execution. The reduced length of an X*-to-X* path Q ^ d[, d2, ■ ■ ■ ,d'j in H* is 



{£(4) - 0[headG. id',)] + c/.[tailG* (d',)]) = ^ £{d[) - 0[end(Q)] + 0[start(Q)]. (6) 



i=l 



Therefore, for any nodes x,y G X* , the x-to-y distance in H* w.r.t. the residual capacity is given by 
D* [x, y] — (f>[y] + (f>[x]. Since the procedure maintains the restriction of <p to faces of X* , this distance can be 
obtained in constant time. Adapting FR's data structure to handle reduced distances w.r.t. a price function 
4> was discussed in Section [3.2. 1[ 



Algorithm 4 Efhcicnt Implementation of hue 11 of AbstractFixConservationOnPath (Algorithm 2) 



% push flow on d to make its residual capacity zero as required for Hassin's algorithm 
1: residual capacity :— c[d] — f[d] — (p [head^. (d)] + (p [tailc. (d)] 

2: val := min{residual capacity, \v[pi\\} amount of flow pushed on d 

3: f[d] = f[d] + val ; /[rev(d)] := -f[d] 

4: w[tail((i)] = 7;[tail(d)] — val ; T;[hcad((i)] = v[head((i)] + val update the inflow at pi and pi+i 

% push excess inflow from tail(d) to head(d) using Hassin's algorithm 

5: let £\d'] ~ C\d'] - f\d'] for all darts d' G PU rev(P) ^^f^'' of explicit darts are residual eapaeities 

L J L J ./ L J \ / (not including circulation component) 

6: £[rev(d)] := \vM - 0^ [tail^. (rev(d))] + 0^ [headc- (rev(d))] r^s.^X" 

„ ( \ FHiY)^ £ (h^ head ffi)) face potential are distances in G* from lieadf^* (d) w.r.t. the 

■ 'Pi \ ) ■— \ : :<? : G* \ ) ) reduced lengths induced by </)^ 

val := 0f [head(d)] - </.f [tail(d)] the amount of flow assigned to d by the circulation corresponding to 

f[d\ := f[d] — val ; /[rev((i)] := —f[{d)] do not push the circulation on d and rev(d) 

accumulate the current circulation 

?;[tail((i)] — i;[tail(c?)] — val ; t;[head((i)] = t;[head((i)] + val update the inflow at pt and pi+i 



It follows from the above discussion that, after executing line H of Algorithm [2] using the efficient imple- 
mentation (Algorithmic]) for the last time, the flow assignment / maintained by the efficient implementation 
is the same as the one that would have been computed if the inefficient implementation (Algorithm [3| were 
used. Moreover, the potential function cp'^ computed by the efficient implementation is the restriction to 
X* of the potential function (p that would have been computed by the inefficient implementation. 

In order to output the flow /' in the entire graph G, we need to know the potential function <p rather 
than just its restriction cj)'^ . Observe, however, that any pseudoflow that differs from / by a circulation is a 
valid output of FixConservationOnPath since a circulation does not change the inflow at any node, nor 
does it introduce residual paths between nodes that are not connected by a residual path in . It therefore 
suffices to find any feasible circulation x iii Gf. This can be done by computing shortest paths in G* w.r.t. 
the lengths £ .= c — f from an arbitrary node x d G* . Note, however, that for darts of P U rev(P) these 
lengths might be negative. We therefore use the 0{n log^ n/ log log n)-time algorithm for shortest paths with 
negative lengths in planar graphs [30] to compute a feasible circulation x- 

The pseudocode for the efficient implementation of FixConservationOnPath is given in appendix [P] 



Running Time Analysis Let n and k be the number of nodes in G and in P, respectively. The initial- 
ization time is dominated by the 0(nlogn -f k^Xogn) time for MSSP. The execution of each iteration of 
the main loop is dominated by the call to FR, which takes 0(A:log^A:) time. The number of iterations is 
k — 1, leading to a total of 0{k^ log^ k) time for execution of the main loop. Computing the circulation x 
requires one shortest path computation, which takes 0(nlog^ n/ log log n) time. Thus total running time of 
the efficient implementation of FixConservationOnPath is 0(nlog^ n/ log log n + fc^ log^ k). 

■''In appendix |e] we show how to use one more call to FR followed by a call to Dijfcstra's algorithm to compute these 
distances. While we think doing so is more elegant and simpler to implement, it does not change the asymptotic running time 
of the algorithm. 
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4 Pushing Excess Inflow from a Cycle 



In this section we describe the procedure CycleToSingleSinkLimitedMaxFlow that pushes excess in- 
flow from a cycle to a node not on the cycle. The procedure is given as Algorithm [5] The procedure 
SingleSourceToCycleLimitedMaxFlow is very similar. We omit its description. 



Algorithm 5 CycleToSingleSinkLimitedMaxFlow(G, c, /o, C, i) 

Input: a directed planar graph G with capacities c, a pseudoflow /o, a simple cycle C, a sink a^. 

Assumes: Vw e C+, infiow^„(u) > 0, where C+ = {v € C : {x ^ C : inflow^, (x) > 0} 4° v}. 
Output: a pseudoflow / s.t. (i) f — fo obeys conservation everywhere but C+ U {t}, 

(a) Vw e C+, inflow/(t;) > 0, (in) {v € C : infiowj^(i;) > 0} Si a,,. 

1: let C+ := {w e C : there exists a residual path to v from a node x € C with inflow (a;) > 0} 
2: delete the nodes of C — C+ 

3: let t;i , W2 , • • • , be the nodes of C+ , labeled according to their cyclic order on C 

4: add artificial arcs WiWi+i for \ < i < t 

5: let P be the vi-to-vi path of artificial arcs 

6: contract all the edges of P, collapsing C+ into the single node Vi 

7; / SiNGLESOURCeSiNGLESiNKMAXFlOw(G', C - /o,Ci,t) 

8: undo the contraction of the edges of P 

9: / := FixConservationOnPath(G', P C, /o + /) -/o 

10; modify / to push back flow to nodes of C+ whose inflow w.r.t /o + / is negative 

11: / /o + / 
12; return / 



To compute C"*" in Line[T] consider the residual graph of G w.r.t. /q. Add a node v and non-zero capacity 
arcs vw for every node w whose inflow w.r.t. /o is positive (these arcs may not preserve planarity). In 0(|G|) 
time, find the set X of nodes that are reachable from v via darts with non-zero capacity. Then C+ ^ C C^X. 

Since C"*" consists of all nodes of C reachable via residual paths from the nodes of C with positive infiow, 
the flow computed by the procedure involves no darts incident to nodes in C — C"*". Thus, restricting the 
computation to the graph obtained from G by deleting the nodes in C — C+ (line [2| does not restrict the 
computed flow. After deletion, adding artificial arcs between consecutive nodes of (line|4| will not violate 
planarity. Contracting the artificial arcs effectively turns into a single node vi. Next, the procedure 
computes a maximum flow / from C+ to t w.r.t. residual capacities c— /q. This is done by invoking a single- 
source single-sink maximum flow algorithm ^ with source vi and sink t (line[7|. Uncontracting the artificial 
arcs turns / into a maximum C^-io-t fiow in G w.r.t. the residual capacities c — /g. However, some of the 
nodes of C+ may have negative inflow w.r.t. fo + f. In line[9j the procedure calls FixConservationOnPath 
to reroute the flow / among the nodes of so that, w.r.t. fo + f, there are no residual paths from nodes 
of C"*" with positive inflow to nodes of C+ with negative infiow. This implies that any node of C+ that still 
has negative inflow has pushed too much flow. Line [TO] modifies / to push back such excess flow so that no 
node of C"'" has negative infiow w.r.t. /q -|- /. This is done using the procedure of Section [O] described in 
the appendix. Finally, the procedure returns fo + f. See Appendix [F| for a formal proof of correctness. 

We next analyze the running time of this procedure on an n-node graph G and a fc-node cycle C. The 
st-maximum flow computation in line[7]takes 0{n\ogn) time using the algorithm of Borradaile and Klein j^J. 
The running time of the procedure is therefore dominated by the call to FixConservationOnPath in line[9] 
which takes 0{n log^ n/ log log n + log^ k) time. 
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A Correctness of Algorithm [T] 



Lemma A.l. Let p he a circulation. Let u and v be nodes in a graph G. Then 

G ^ 

Proof. Pushing a circulation does not change the amount of flow crossing any cut. This implies that if there 
was no u-to-v residual path before p was pushed, then there is none after p is pushed. □ 

We will use the following two lemmas in the proof of correctness: 

Lemma A. 2. (sources lemma) Let f be a flow with source set X. Let A, B be two disjoint sets of nodes. 
Then 

AUX A^ B. 

Proof. The flow / may be decomposed into a cyclic component (a circulation) and an acyclic component. 



By Lemma A.l it suffices to show the lemma for an acyclic flow /. 

Suppose for the sake of contradiction that there exists a residual a-to-& simple path P after / is pushed 
for some a E A and b E B. Let P' be the maximal suffix of P that was residual before the push. Maximality 
implies that the dart d of P whose head is start(P') was non-residual before the push, and f{rev{d)) > 0. 
The fact that f{rev{d)) > implies that before / was pushed there was a residual path Q from some node 
X £ X to head((i). Therefore, the concatenation of Q and P' was a residual a;-to-6 residual path before the 
push, a contradiction. □ 

The proof of the following lemma is symmetric. 

Lemma A. 3. (sinks lemma) Let f be a flow with sink set X . Let A, B be two disjoint sets of nodes. Then 

A^B\JX^ A^B. 

For node sets Y, Z, we will use the notation sources lemma{W, Y, Z) to indicate the use of the sources 
lemma to establish that there are no W-to-Y residual paths after a flow with source set Z is pushed. We 
will use sinks lemma{W,Y, Z) in a similar manner. 

Lemma A. 4. At any time in the running of the algorithm after the last execution of line and before the 
execution of line \13\ S~^T,S^A, S^CjA^TjC-y^T. 

Proof. The fact that the properties hold just after line [7] follows from the properties of the recursive calls and 
from the fact that any residual path from a node in one piece to a node in the other consists of a residual 
path to C and a residual path from C. Note that, because of the contractions in line[4j at this time the cycle 
C consists of just the node v. The the nonexistence of residual paths with respect to v in the recursive calls 
implies the nonexistence of residual paths w.r.t. any node of C after the contractions are undone in linejSj 



Now, let us prove that the properties are maintained until just before the execution of line 13 By the 
above argument, there is a cut separating S from TUAUC which is saturated just after linejSj The procedure 
in line |9] only pushes flow between vertices of C, and in lines \T2\ flow is only pushed between the nodes 
of A and C. These sets are all on the same side of the cut which therefore stays saturated. It follows that 
S ^ T, S ^ A, and S" C at any point after line [7] and before line 13 A similar argument applied to a 



saturated cut separating A(J C from T shows A T and C -f* T. □ 

Recall that Cf' is the set of nodes of C that are reachable via residual paths from some node of C with 
positive excess at the beginning of iteration i, and that C~ is the set of nodes of C that have residual paths 
to some node of C with negative excess at the beginning of iteration i. 

Lemma A. 5. Just after /me[^ of Algorithm^^is executed, Cf . 

Proof. Follows from the definition of FixConservationOnPath. □ 
Lemma A.6. For all i < j, C+ C C+ and Cj C C~ 
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Proof. It suffices to show that, for aU i, Cj,^ C and C~j^^ C C~ . 

The flow pushed in iteration i of hue 11 can be decomposed into a flow whose sources and sinks are aU 
in and a flow whose sources are in and whose sink is a^. Therefore, the set X of nodes of C with 
positive inflow immediately after iteration i of line 11 is a subset of Cf . By definition of SingleSource- 



ToCycleLimitedMaxFlow, the set of nodes of C with positive inflow does not change after hue [12| is 
executed. Therefore, Cj^^ is the set of nodes reachable from X by a residual path after iteration i of line|l2l 
By definition of , immediately before iteration i of line 11 there are no C,^-to-{C — C^} residual 



paths. By sources lemma{Cf ^ C — Cf , C^), there are no Cj^-to-{C — C^} residual paths immediately after 
iteration i of line 11 as well. This shows that there are no X-to-{C — C^} residual paths at that time. 

can be decomposed into a flow whose sources and sinks are all in and 



The flow pushed in line 
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a flow whose source is and whose sinks are all in . By sinks le mm a{Gf ^ C — Cf ), there are 

This shows that there are no 
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no Cf-to-{C — C^} residual paths immediately after iteration i of line 
X-io-{C — cf} residual paths at that time. Hence Cf^^ C Cf , as desired. 
The proof of the analogous claim for Cj^j^ is similar. 



□ 



Lemma A. 7. Just after line 11 of Algorithm^^ is executed in iteration i, Cf -f* , Cf^^ -f^ \fij\i<iy 



j )j<i 



Proof. The flow pushed in line 11 can be decomposed into a flow whose sources and sinks are all in Cf and 
a flow whose sources are in Cf and whose sink is a^. 

• cf C^ by sources lemma{Cf,C^,Cf ) 

• ^i+i ^ {%}j<i by sources lemma{Cf^i,aj,Cf) 



C 



i+l 



{oi} by definition of CycleToSingleSinkLimitedMaxFlow 



{aj}j<i ^ Cj by sources lemma{aj,C^ ,Cf ) 



Lemma A. 8. Just after line 12 of Algorithm [7] is executed in iteration i, C^ 



□ 



{aj}]<i 



C7, C 



i+i 



j Sj<i! 



C 



■i+l- 



Proof. The flow pushed in line 12 can be decomposed into a flow whose sources and sinks are all in C, and 



a flow whose source is and whose sinks are all in C~ . 

• cf C^ by sinks leninia{Cf , C~, C^) 

• ^i+i ^ {%}j<'i t>y sinks lemma{Cfj^^,aj^C~) 

• {%'}j<'i ^ ^i+i sinks lemma{aj,C^_^-^,C^) 

• ai ^ Cj^j by definition of SingleSourceToCycleLimitedMaxFlow 



Let C+(C^) denote the set of nodes in C with positive (negative) inflow just before line 13 
Corollary A. 9. Just before line\l^ of Algorithm\]\ is executed, C+ C^,C^ A,A^ C^ . 
Lemma A. 10. The following are true upon termination: 

1. f is a pseudoflow 

2. f obeys conservation everywhere except at S, T, A 

3. s^T,s'^ A,A^T. 



□ 

is executed. 
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Proof. Since every addition to / along the algorithm respects capacities of all darts, / is a pseudoflow at 
all times. By induction, the only nodes that do not obey conservations after the recursive calls are those of 
S, T and A. Subsequent changes to / only violate conservation on the nodes of C, but any such violation is 
eliminated in line |13| Therefore upon termination / obeys conservation everywhere except S, T and A. 



Since, by Lemma A. 4 and Corollary A. 9 before line 13 C+ ^ A and C ^ T, the flow pushed back from 
in line 13 must be pushed back to S. Similarly, the flow pushed back to C~ must be pushed back from 

T. 



Let /+ (/_) be the flow pushed back from C+ to S (from T to C ) in line 13 Consider first pushing 
back /+. 

• 5 T by sources lemma{S, T, C^) 

• S ^ Ahy sources lemma{S, A, C^) 
A T hy sources lemma{A, T, C+) 

m S C~ by sources lemma{S, C~, C^) 
A by sources lemma{A, C^, C'^) 

Next consider pushing /_ 

• 5 T by sinks lemma{S, T, C~) 

• 5 ^ j4 by sinks lemma{S, A, C^) 

• A-^ T \>y sinks lemma{A, T, C") 

□ 

B Correctness of Algorithm [2] 

Lemma B.l. The following holds immediately after iteration i of the main loop (line^. 

1. For j < i,j' > j, if Pj has positive inflow, there is no residual path from pj to pji . If Pj has negative 
inflow, there is no residual path from pj> to pj . 

2. For < i, if Pj has positive inflow and pj' has negative inflow then there is no pj-to-pj> residual 
path. 

Proof. By induction on the number of iterations i of the loop. For i = the invariants are trivially satisfied. 

First note that the adjustments to capacities and fiow in lines |6j|8] do not create any new residual paths 
since capacities are only reduced, and no residual capacity increases. Therefore, it suffices to argue just 
about the flow pushed in line[TT| 

Assume the invariants hold up until the beginning of the i*'* iteration. We show that the invariants hold 
at the end of the iteration. Suppose that pi has positive inflow at the end of the i*^ iteration (the case of 
negative inflow is similar). 



1. Since the flow pushed in line 11 is limited by |inflow(pj)|, the fact that pi has positive inflow at the 
end implies that the flow pushed was a maximum flow from pi to Pi+i. Since the capacities of darts 
between pk+i and pk for k > i are sufficiently large, the maximality of the fiow implies that there are 
no pi-to-pk residual paths for any k > i. 

The invariant holds for nodes pj with positive infiow and j < i hy sinks lemma{{pj}, {pji : j' > j}, 

te+i})- 

2. Invariant [2] holds for < ? by sinks lemma{ {pj : j < i,m&ow{pj) > 0}, {pj : j < i,inflow(pj) < 0}, 

te+i})- 

The invariant holds for pi by invoking the sources lemma{{pi}, {pj : j < i,inflow(pj) < 0}, {pi})- 
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□ 



Corollary B.2. Algorithm^ is correct 

Proof. The flow /' — /o satisfies conservation everywhere except at nodes of P since the algorithm only 
pushes flow between nodes of P. By Lemma B.l with respect to /', there are no residual paths from nodes 
of P with positive inflow to nodes of P with negative inflow. □ 



C Converting a Pseudoflow into a Maximum Feasible Flow 

Let / be a pseudoflow in a planar graph G with node set V, sources S and sinks T. Let denote the set 
of nodes {v eV ~ (SUT) : inflow/(u) > 0}. Similarly let V' denote the set of nodes {v eV - (SUT) : 

inflow J (f) < 0}. Suppose S U T U . In this appendix we show how to convert / into a maximum 

feasible flow /'. This procedure was first described for planar graphs by Johnson and Venkatesan ^9J. The 
original description of the procedure took O(nlogri) time, but using the flow cycles canceling technique of 
Kaplan and Nussbaum [20j the running time is 0{n). 

We begin by converting / to an acycic pseudoflow in linear time [20J. That is, after the conversion there 
is no cycle C such that f[d] > for every dart d of C. Since / is acyclic, the darts with a positive flow 
induce a topological ordering on the nodes of the graph G. Let v be the last member of in the topological 
ordering, and let d be an arbitrary dart that carries flow into v. We set f[d] — max{/[(i] — inflowj^(w), 0}, 

and set /[rev(d)] accordingly. The flow assignment / maintains the invariant S U T U V~ by sinks 

lemma{S U ,T U V~ ,{v}). As long as v is in there must be a dart d which carries flow to v. By 
changing the flow on d we cannot add to V'^ a new node that appears later than v in the topological ordering. 
We repeat this process until is empty. Since we reduce the flow on each dart at most once, this takes 

linear time. Next we handle V~ while keeping the invariant S U -d T U in a symmetric way, by 
repeatedly fixing the first vertex if in the topological ordering. 

The total running time is 0(n), and since , are both empty, we get from the invariant S U 
T U that the resulting pseudoflow is a feasible flow. This is the required flow /'. 
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D Pseudocode for FixConservationOnPath 



Algorithm 6 FixConservationOnPath(G, P, c, /o) 



Input: a directed planar graph G, simple path P = did2 ■ ■ - dk, capacity function c, and a pseudoflow /o 
Output: a pseudoflow / s.t. (1) f — fo satisfies conservation at nodes not on P and (2) w.r.t /, no residual 
path exists from {v € P : inflow(i;) > 0} to {u G P : inflow(w) < 0}. 

let X* be the set of endpoints in the planar dual G* of the darts of P 



initialize to zero an array <p'^ [■] indexed by the nodes of X* 
let H* be the graph obtained from G* by removing the darts of P 

compute, using the MSSP algorithm [23l[5], a table D*[-, •] of distances in H* between nodes oi X* w.r.t 
the lengths c — fo 
/:=/o 

c[d] := c[d] + M for all darts of P U rev(P) 

initialize an array of length k by v[i] := inflow at pi with respect to / 
for i — 1,2, . . . , k 

let Pi and Pi+i be the tail and head of di, respectively 

% reduce the capacities of d and rev(d) by M and adjust the flow appropriately 
10: for d e {di,rev{di)} 
11: c[d] := c[d] - M 

12: old flow := f[d] + ^^[headG*(d)] - 0^[tailG* (d)] 

13: new flow := min{old flow, c[d]} adjusted flow must not exeeed adjusted capacity 

14: f[d] :— new flow — 0^[headG» {d)] + (j)'^ [tailQ* (d)] the exphcit flow does not include the circulation component 

15: /[rev(d)] -/[d] 

16: t>[head(o?)] :— v[head{d)] + new flow — old flow update the inflow at pt and 

17: t>[tail((i)] := t;[tail(d)] — new flow + old flow 

18: if v[pi\ > then d := di else d := Tev{di) flnd in which direction flow should be pushed 

% push flow on d to make its residual capacity zero as required for Hassin's algorithm 
19: residual capacity := c[d\ — f[d] — 0^ [headc* (d)] + 0'^[tailG. (d)] 

20: val := min{residual capacity, \v[pi]\} amount of flow pushed on d 

21: m = f[d] + val ; /[rev(d)] -f[d] 

22: t)[tail(d)] = t>[tail(d)] — val ; i;[head(d)] = D[head(d)] + val update the inflow at p; and pi+i 

% push excess inflow from tail(d) to head(ci) using Hassin's algorithm 

23: let i\d'] Cfd'l - fid'] for aU darts d' € P U rev(P) of explicit darts are residual capacities 

L -I L J ./ L J \ / (not mcludmg circulation component) 

24: 4rev(d)]:=mp.]|-0^[tailG.(rev(d))] + 0^[headG*(rev(d))] rad^u^X"™^^^ 

25" cb^ ( ) FR(Y)^ £ d)^ head (d)') face potential are distances in G* from header* (d) w.r.t. the 

''■ ^ ' ' ^ 1 •> ^ 1 \ // reduced lengths induced by cf>'^ 

val := (t}f [head(rev(d))] - </>f [tail(rev(d))] the amount of flow assigned to d by the circulation corresponding to (p^ 



26: 
27: 
28: 
29: 

30: 
31: 
32: 
33: 
34: 



f[d\ :— f[d\ + val ; /[rev(d)] :— —/[(d)] do not push the circulation on d and rcv(ci) 

accumulate the current circulation 

t;[tail(d)] = t;[tail(d)] — val ; t;[head(d)] — t)[head(d)] + val update the inflow at pi and pi+i 

% find a feasible circulation 

let X be an arbitrary node in G* 

let £[d] := c[d] - f[d] for all darts d € G 

X ■= SingleSourceShortestPaths(G*,^, x) distances of the nodes of G* from x 

/'[d] := f[d\ + x[head(3* (d)] — x[tailG'» (d)] for every dart d € G add explicit flow and implicit circulation 
return /' 
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E Alternative to Line [32] of Algorithm [6] 



In this section we show that it is not necessary to use a generic shortest path algorithm that works with 
negative lengths to compute a feasible circulation x in line [32] of Algorithm |6] Instead of choosing x to be 
an arbitrary node in G* , let x be an arbitrary node of X* . Let xiv) denote the x-to-y distance in G* w.r.t. 
the lengths £. Instead of computing x using a shortest path algorithm that accepts negative lengths, we will 
compute it more efficiently in two steps. Pseudocode is given below as Algorithm [7j In the first step we use 
FR to compute the distances to just the nodes of X* . In the second step we extend these distances to all 
other nodes using Dijkstra's algorithm. 

In the first step the algorithm computes distances in G* from x to all nodes of X* w.r.t. £^x , the reduced 
lengths of £ induced by the feasible price function (p'^ . This is done by an additional invocation of FR 
(line[2|. Let xp^ denote these distance labels. By definition of reduced lengths and a telescoping sum similar 
to the one in the derivation of Eq. ^ , 

xp''[y]^x{y) + <t>''[^]-cp''[y]- (7) 

Since both (p'^ and tp'^ are known, the algorithm can compute the unreduced distances xiu] foi' U G ^* 
(line|4|. Next, the algorithm runs Dijkstra's algorithm on G* , initializing the label of each node y of X* to 
its correct value xiu]- Since in the dual the darts of P U rev(P) are only incident to nodes of X* , Dijkstra's 
algorithm initialized in this manner correctly outputs the distance labels for all nodes of G* even if some of 
the darts of P U rev(P) may have negative lengths. 



Algorithm 7 Alternative to line 32 of Algorithm |6] 



let x be an arbitrary node in A"* 

find distances of tlic nodes of X* from x w.r.t. price function <p 

initialize to oo an array x indexed by the nodes of G* 

X[y] = '0'^ [y] ~ 'P^ l"^] + 4>^ [y] every y € X* distances of tlie nodes of X' from x 

extend £ to all darts of G* by £[<]] c[d\ — fo[d] for every dart d £ G 

X ■= Dijkstra(G'*,^, X, distances of the nodes of G* from x 



F Correctness of Algorithm [5] 

Lemma F.l. Algorithm^is correct. 

Proof. Any flow that is pushed by the algorithm originates at and terminates at C+ U {t}. Therefore, 



/ — /o violates conservation only at U {t}. In line 10 flow of / is pushed back so that no node of C+ 
has negative inflow w.r.t. /o + /. It is possible to do so by only pushing back flow of / (rather than flow of 
fo + /) since by assumption no node of has negative inflow w.r.t. /q. 

By maximality of the flow pushed in line [?) just after line [7| is executed there are no C+-to-i residual 
paths. Clearly, this remains true when the capacities of the artificial darts are set to zero. In line [9] flow is 
pushed among the nodes of C"*", so by .sources lemma(C^ ,t,C^), there are no C"'"-to-i residual paths after 
line|9]either. Moreover, by definition of FixConservationOnPath, there are no C>o-to-C<o residual paths 
immediately after line [9] is executed, where C>o (C'<o) is the set of nodes of C+ with positive (negative) 
inflow at that time. Line |10| pushes flow into C<o, making all the nodes of C<o obey conservation. By sinks 
lemma{CyQ, t, C<o) there are no C^Q-to-t residual paths upon termination of the procedure. This completes 
the proof since C>o is the set of nodes of C with positive inflow upon termination. □ 
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